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Abstract. We introduce two simple models of forward-backward stochastic differential equations 
with a singular terminal condition and we explain how and why they appear naturally as models for the 
valuation of CO2 emission allowances. Single phase cap-and-trade schemes lead readily to terminal 
conditions given by indicator functions of the forward component, and using fine partial differential 
equations estimates, we show that the existence theory of these equations, as well as the properties of 
the candidates for solution, depend strongly upon the characteristics of the forward dynamics. Finally, 
we give a first order Taylor expansion and show how to numerically calibrate some of these models for 
the purpose of CO2 option pricing. 



1. Introduction 

This paper is motivated by the mathematical analysis of the emissions markets, as implemented 
for example in the European Union (EU) Emissions Trading Scheme (ETS). These market mecha- 
nisms have been hailed by some as the most cost efficient way to control Green House Gas (GHG) 
emissions. They have been criticized by others for being a tax in disguise and adding to the burden 
of industries covered by the regulation. Implementation of cap-and-trade schemes is not Umited to 
the implementation of the Kyoto protocol. The successful US acid rain program is a case in point. 
However, a widespread lack of understanding of their properties, and misinformation campaigns by 
advocacy groups more interested in pushing their political agendas than using the results of objective 
scientific studies have muddied the water and add to the confusion. More mathematical studies are 
needed to increase the understanding of these market mechanisms and raise the level of awareness of 
their advantages as well as their shortcomings. This paper was prepared in this spirit. 

In a first part, we introduce simple single-firm models inspired by the workings of the electricity 
markets (electric power generation is responsible for most of the CO2 emissions worldwide). Despite 
the specificity of some assumptions, our treatment is quite general in the sense that individual risk 
averse power producers choose their own utility functions. Moreover, the financial markets in which 
they trade emission allowances are not assumed to be complete. 

While market incompleteness prevents us from identifying the optimal trading strategy of each 
producer, we show that, independently of the choice of the utility function, the optimal production or 
abatement strategy is what we expect by proving mathematically, and in full generality (i.e. without 
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assuming completeness of the markets), a/o/fc theorem in environmental economics: the equilibrium 
allowance price equals the marginal abatement cost, and market participants implement all the abate- 
ment measures whose costs are not greater than the cost of compliance (i.e. the equilibrium price of 
an allowance). 

The next section puts together the economic activities of a large number of producers and searches 
for the existence of an equilibrium price for the emissions allowances. Such a problem leads naturally 
to ?i forward stochastic differential equation (SDE) for the aggregate emissions in the economy, and 
a backward stochastic differential equation (BSDE) for the allowance price. However, these equa- 
tions are "coupled" since a nonlinear function of the price of carbon (i.e. the price of an emission 
allowance) appears in the forward equation giving the dynamics of the aggregate emissions. This 
feedback of the emission price in the dynamics of the emissions is quite natural. For the purpose of 
option pricing, this approach was described in [5] where it was called detailed risk neutral approach. 

Forward backward stochastic differential equations (FBSDEs) of the type considered in this section 
have been studied for a long time. See for example |[T4l . or |[T9l . However, the FBSDEs we need 
to consider for the purpose of emission prices have an unusual pecularity: the terminal condition 
of the backward equation is given by a discontinuous function of the terminal value of the state 
driven by the forward equation. We use our first model to prove that this lack of continuity is not 
an issue when the forward dynamics are strongly elliptic in the neighborhood of the singularities 
of the terminal condition, in other words when the volatility of the forward SDE is bounded from 
below in the neighborhood of the discontinuities of the terminal value. However, using our second 
equilibrium model, we also show that when the forward dynamics are degenerate (even if they are 
hypoelliptic), discontinuities in the terminal condition and lack of strong ellipticity in the forward 
dynamics can conspire to produce point masses in the terminal distribution of the forward component, 
at the locations of the discontinuities. This implies that the terminal value of the backward component 
is not given by a deterministic function of the forward component, for the forward scenarios ending 
at the locations of jumps in the terminal condition, and justifies relaxing the definition of a solution 
of the FBSDE. 

Even though we only present a detailed proof for a very specific model for the sake of definiteness, 
we believe that our result is representative of a large class of models. Since from the point of view 
of the definition of "aggregate emissions", the degeneracy of the forward dynamics is expected, this 
seemingly pathological result should not be overlooked. Indeed, it sheds new light on an absolute 
continuity assumption made repeatedly in equilibrium analyses, even in discrete time models. See for 
example f45 and fi \. This assumption was regarded as an annoying technicality, but in the light of the 
results of this paper, it looks more intrinsic to these types of models. In any case, it fully justifies the 
need to relax the definition of a solution of a FBSDE when the terminal condition of the backward 
part jumps. 

A vibrant market for options written on allowance futures/forward contracts has recently developed 
and increased in liquidity. See for example |5| for details on these markets. Reduced formed models 
have been proposed to price these options. See |5| or 161. Several attempts have been made at 
matching the smile (or lack thereof) contained in the quotes published daily by the exchanges. Section 
[5] develops the technology needed to price these options in the context of the equilibrium framework 
developed in the present paper. We identify the option prices in terms of solutions of nonlinear partial 
differential equations and we prove when the dynamics of the aggregate emissions are given by a 
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geometric Brownian motion, a Taylor expansion formula when the nonlinear abatement feedback is 
small. We derive an explicit integral form for the first order Taylor expansion coefficient which can 
easily be computed by Monte Carlo methods. We believe that the present paper is the first rigorous 
attempt to include the nonlinear feedback term in the dynamics of aggregate emissions for the purpose 
of emissions option pricing. 

The final Section|5]was motivated by the desire to provide practical tools for the efficient computa- 
tion of option prices within the equilibrium framework of the paper. Indeed, because of the nonlinear 
feedback created by the coupling in the FBSDE, option prices computed from our equilibrium model 
differ from the linear prices computed in 0, fT\\ and HI in the framework of reduced form models. 
We derive rigorously an approximation based on the first order asymptotics in the nonlinear feedback. 
This approximation can be used to compute numerically option prices and has the potential to effi- 



ciently fit the implied volatility smile present in recent option price quotes. The final subsection 5.3 
illustrates numerically the properties of our approximation. 

Acknowledgements: We would like to thank two anonymous referees for pointing out inconsisten- 
cies in the original proofs of Theorem [T] and Proposition [4| and for insightful comments which led to 
improvements in the presentation of the results of the paper. 



2. Two Simple Models of Green House Gas Emission Control 

We first describe the optimization problem of a single power producer facing a carbon cap-and- 
trade regulation. We assume that this producer is a small player in the market in the sense that his 
actions have no impact on prices and that a liquid market for pollution permits exists. In particular, we 
assume that the price of an allowance is given exogenously, and we use the notation Y = (lt)o<t<T 
for the (stochastic) time evolution of the price of such an emission allowance. For the sake of sim- 
plicity we assume that [0, T] is a single phase of the regulation and that no banking or borrowing of 
the certificates is possible at the end of the phase. For illustration purposes, we analyze two simple 
models. Strangely enough, the first steps of these analyses, namely the identifications of the optimal 
abatement and production strategies, do not require the full force of the sophisticated techniques of 
optimal stochastic control. 

2.1. Modeling First tlie Emissions Dynamics. We assume that the source of randomness in the 
model is given hyW = {Wt)Q<t<T, a finite family of independent one-dimensional Wiener processes 
= {Wl)o<t<T, l<3<d.\n other words, Wt = iW^ , ■■■ , Wf) for each fixed t G [0, T]. All 
these Wiener processes are assumed to be defined on a complete probability space J^, F), and we 
denote by F = {J-f , t > 0} the Brownian filtration they generate. Here, T > is a fixed time horizon 
representing the end of the regulation period. 

We will eventually extend the model to include N firms, but for the time being, we consider only 
the problem of one single firm whose production of electricity generates emissions of carbon dioxyde, 
and we denote by Et the cumulative emissions up to time t of the firm. We also denote by Et the 
perception at time t (for example the conditional expectation) of what the total cumulative emission 
Et will be at the end of the time horizon. Clearly, E and E can be different stochastic processes, but 
they have the same terminal values at time T, i.e. Et = Et- We will assume that the dynamics of 
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the proxy E for the cumulative emissions of the firm are given by an Ito process of the form: 

(1) Et = Eo+ [ {bs - ^s)ds + [ asdWs, 

Jo Jo 

where b represents the (conditional) expectation of what the rate of emission would be in a world 
without carbon regulation, in other words in what is usually called Business As Usual (BAU for 
short), while ^ is the instantaneous rate of abatement chosen by the firm. In mathematical terms, ^ 
represents the control on emission reduction implemented by the firm. Clearly, in such a model, the 
firm only acts on the drift of its perceived emissions. For the sake of simplicity we assume that the 
processes b and a are adapted and bounded. Because of the vector nature of the Brownian motion W, 
the volatility process a is in fact a vector of scalar volatility processes {cr^)i<j<d- For the purpose 
of this section, we could use one single scalar Wiener process and one single scalar volatility process 
as long as we allow the filtration F to be larger than the filtration generated by this single Wiener 
process. This fact will be needed when we study a model with more than one firm. 

Notice that the formulation ([T]) does not guarantee the positiveness of the perceived emissions 
process, as one would expect it to be. This issue will be discussed in Proposition [3] below, where we 
provide sufficient conditions on the coefficients of ([T]l in order to guarantee the positiveness of the 
process E. 

Continuing on with the description of the model, we assume that the abatement decision is based 
on a cost function c : M — M which is assumed to be continuously differentiable (C^ in notation), 
strictly convex and satisfy Inada-like conditions: 

(2) c'(— oo) = — oo and c'(+oo) = +oo. 

Note that (c')~^ exists because of the assumption of strict convexity. Since c{x) can be interpreted 
as the cost to the firm for an abatement rate of level x, without any loss of generality we will also 
assume c(0) = min c = 0. Notice that ^ implies that lim^^-^ioo c{x) = +oo. 

Example 1. A typical example of abatement cost function is given by the quadratic cost function 
c{x) = ax"^ for some a > used in [21], or more generally the power cost function c{x) = 
for some a > and /3 > 0. 

The firm controls its destiny by choosing its own abatement schedule ^ as well as the quantity 6 of 
pollution permits it holds through trading in the allowance market. For these controls to be admissible, 
^ and 6 need only to be progressively measurable processes satisfying the integrability condition 

(3) E [ [e^ + ^f]dt < oo. 

Jo 

We denote by A the set of admissible controls 9). Given its initial wealth x, the terminal wealth 
Xt of the firm is given by: 

(4) Xt = = X + [ OtdYt - [ c{^t)dt - EtYt. 

Jo Jo 

The first integral in the right hand side of the above equation gives the proceeds from trading in the 
allowance market. Recall that we use the notation Yt for the price of an emission allowance at time 
t. The next term represents the abatement costs, and the last term gives the costs of the emission 
regulation. Recall also that at this stage, we are not interested in the existence or the formation of this 
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price. We merely assume the existence of a liquid and frictionless market for emission allowances, 
and that Yt is the price at which each firm can buy or sell one allowance at time t. The risk preferences 
of the firm are given by a utility function [/ : M — M, which is assumed to be C^, increasing, strictly 
concave and satisfying the Inada conditions: 

(5) {Uy {-oo) = +00 and (?7)'(+oo) = 0. 
The optimization problem of the firm can be written as the computation of: 

(6) V{x) = sup EU{X^^) 

where E denotes the expectation under the historical measure P, and A is the set of abatement and 
trading strategies 6) admissible to the firm. The following simple result holds. 

Proposition 1. The optimal abatement strategy of the firm is given by: 

Ct = [cTHYt). 

Remark 1. Notice that the optimal abatement schedule is independent of the utility function. The 
beauty of this simple result is its powerful intuitive meaning: given a price Yt for an emission al- 
lowance, the firm implements all the abatement measures which make sense economically, namely all 
those costing less than the current market price of one allowance (i.e. one unit of emission). 

Proof. By an immediate integration by parts in the expression Q of the terminal wealth, we see that: 

EtYt = Yt (eo + ^ hdt + ^ atdWtj -^t Ctdt 

= Yt (^Eo + ^ hdt + ^ atdW?j - ^ Yt^tdt - (^j^ ^sds^ dYt 
so that Xt = At + Bi with 



A^= f OtdYt - Yt (eo + I hdt + [ utdWt 
Jo \ Jo Jo 



where the modified control 9 is defined by dt = 9t + /q £,sds, and 



4 = x- r[c{it)-Ytit]dt. 
Jo 



Notice that depends only upon ^ without depending upon 6 while depends only upon 6 without 
depending upon The set A of admissible controls is equivalently described by varying the couples 
{9, ^) or (6*, ^), so when computing the maximum 

sup m{XT) = sup m{A'^T + B^) 

one can perform the optimizations over 9 and ^ separately, for example by fixing 9 and optimizing 
with respect to ^ before maximizing the result with respect to 9. The proof is complete once we notice 
that U is increasing and that for each t G [0, T] and each uj €i Q,, the quantity B^ is maximized by 
the choice = (c')"H^t)- ^ 
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Remark 2. The above result argues neither existence nor uniqueness of an optimal admissible set 
of controls. In the context of a complete market, once the optimal rate of abatement ^* 
is implemented, the optimal investment strategy 9* should hedge the financial risk created by the 
implementation of the abatement strategy. This fact can be proved using the classical tools of portfolio 
optimization in the case of complete market models. Indeed, if we introduce the convex dual UofU 
defined by: 

U{y) := swp{U{x) - xy} 



X 

T'\-l 



and the function I by I = [U') so that U{y) = U o I(y) — yl{y) and if we denote by E and 'EP 
respectively the expectations with respect to P and the unique equivalent measure Q under which Y 
is a martingale (we write Ztfor its volatility given by the martingale representation theorem), then 
from the a.s. inequality 

dQ .0 rf dQ\ dQ f dQ 

valid for any admissible 6), and y we get 

EU(XI:») <EUol(y'§y yE" [^■» - / (y'§ 

after taking expectations under P. Computing E^^X^^ by integration by parts we get: 

EC/(x|,'^)<EC/o/(^y^j+y x-W!^ j\co{c^)-\Yt) + Yt{ht- {c^r\Yt)) + atZt]dt 

-m (y- 



if we use the optimal rate of abatement. So if we choose y = y G M a* the unique solution of: 

y^) =x-E'^j\o ic')-\Yt) + Ytih - {c')-\Yt)) + atZtdt. 



it follows that 



E«X|'^ = E'O/ ( y 



dF 



and finally, if the market is complete, the claim I (y^^ is attainable by a certain 9*. This completes 
the proof. 

2.2. Modeling the Electricity Price First. We consider a second model for which again, part of the 
global stochastic optimization problem reduces to a mere path-by-path optimization. As before, the 
model is simplistic, especially in the case of a single firm in a regulatory environment with a liquid 
frictionless market for emission allowances. However, this model will become very informative later 
on when we consider N' firms interacting on the same market, and we try to construct the allowance 
price Yt by solving a Forward-Backward Stochastic Differential Equation (FBSDE). The model con- 
cerns an economy with one production good (say electricity) whose production is the source of a 
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negative externality (say GHG emissions). Its price {Pt)o<t<T evolves according to the following Ito 
stochastic differential equation: 

(7) dPt = ti{Pt)dt + a{Pt)dWt 

where the deterministic functions fi and a are assumed to be with bounded derivatives. At each 
time t S [0, T], the firm chooses its instantaneous rate of production qt and its production costs are 
c{qt) where c is a function c : M+ ^ M which is assumed to be and strictly convex. With these 
notations, the profits and losses from the production at the end of the period [0, T], are given by the 
integral: 

[ [Ptqt - c{qt)]dt. 
Jo 

The emission regulation mandates that at the end of the period [0,T], the cumulative emissions of 
each firm be measured, and that one emission permit be redeemed per unit of emission. As before, 
we denote by (lt)o<t<T the process giving the price of one emission allowance. For the sake of sim- 
plicity, we assume that the cumulative emissions Et up to time t are proportional to the production in 
the sense that Et = eQt where the positive number e represents the rate of emission of the production 
technology used by the firm, and Qt denotes the cumulative production up to and including time t: 




At the end of the time horizon, the cost incurred by the firm because of the regulation is given by 
EtYt = eQrYT- The firm may purchase allowances: we denote by 6t the amount of allowances 
held by the firm at time t. Under these conditions, the terminal wealth of the firm is given by: 

(8) Xt = xf = x+ [ OtdYt + / [Ptqt - c{qt)]dt - eQrYT 

Jo Jo 

where as before, we used the notation x for the initial wealth of the firm. The first integral in the 
right hand side of the above equation gives the proceeds from trading in the allowance market, the 
next term gives the profits from the production and the sale of electricity, and the last term gives the 
costs of the emission regulation. We assume that the risk preferences of the firm are given by a utility 
function [/ : M — M, which is assumed to be C^, increasing, strictly concave and satisfying the Inada 
conditions (|5]) stated earlier. As before, the optimization problem of the firm can be written as: 

(9) V{x) = sup EU{xf) 

where E denotes the expectation under the historical measure P, and A is the set of admissible 
production and trading strategies {q,9). This problem is similar to those studied in |[2l where the 
equilibrium issue is not addressed. As before, for these controls to be admissible, q and need only 
be adapted processes satisfying the integrability condition 

fT 

(10) E / [e'f +qf]dt < oo. 

Jo 

Proposition 2. The optimal production strategy of the firm is given by: 

q*t={crHPt-eYt). 



s 
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Remark 3. As before, the optimal production strategy q* is independent of the risk aversion (i.e. the 
utility function) of the firm. The intuitive interpretation of this result is clear: once a firm observes 
both prices Pf and Yf, it computes the price for which it can sell the good minus the price it will have 
to pay because of the emission regulation, and the firm uses this corrected price to choose its optimal 
rate of production in the usual way. 

Proof. A simple integration by part (notice that Et is of bounded variations) gives: 

T 7^ 7^ 

(11) QtYt= [ YtdQt+ [ QtdYt = [ Ytqtdt+ [ QtdYt, 

so that Xt = Ai + Bl with 



^0 ^0 ^0 



A'^=retdYt with ~et = Ot - e I q.^ds 







JO 



which depends only upon 9 and 



B't = x+ [ [{Pt - eYt)qt - c{qt)]dt, 
Jo 



which depends only upon q without depending upon 6. Since the set A of admissible controls is 
equivalently described by varying the couples (q, 9) or {q, 9), when computing the maximum 

sup ¥.{U{Xt)]= sup ¥.{U{A^ + B'^)] 

one can perform the optimizations over q and 9 separately, for example by fixing 9 and optimizing 
with respect to q before maximizing the result with respect to 9. The proof is complete once we notice 
that U is increasing and that for each t G [0, T] and each u £ Q,, the quantity B^ is maximized by 
the choice q^ = (c')"^^* - ^^t). □ 

3. Allowance Equilibrium Price and a First Singular FBSDE 

The goal of this section is to extend the first model introduced in section |2] to an economy with N 
firms, and solve for the allowance price. 

3.1. Switching to a Risk Neutral Framework. As before, we assume that Y = (lj)te[o,T] is the 

price of one allowance in a one-compliance period cap-and-trade model, and that the market for 
allowances is frictionless and liquid. In the absence of arbitrage, y is a martingale for a measure Q 
equivalent to the historical measure P. Because we are in a Brownian filtration. 



exp 



T 1 /"^ ■ 

atdWt — - \at\'^dt 
2 ./n 



for some sequence a = (at)tg[o,T] of adapted processes. By Girsanov's theorem, the process W 

{Wt)telo,T] defined by 



Wt = Wt- I a,ds 
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is a Wiener process for Q so that equation ([T]l giving the dynamics of the perceived emissions of a 
firm now reads: 

dEt = (bt - it)dt + atdWt 
under Q, where the new drift b is defined by ht = bt + crtat for all t G [0, T]. 

3.2. Market Model with Firms. We now consider an economy comprising N firms labelled by 
{!,••• , N}, and we work in the risk neutral framework for allowance trading discussed above. When 
a specific quantity such as cost function, utility, cumulative emission, trading strategy, . . . depends 
upon a firm, we use a superscript * to emphasize the dependence upon the z-th firm. So in equilibrium 
(i.e. whenever each firm implements its optimal abatement strategy), for each firm i E {1, • " " i ^} 
we have 

dEl = m-[{^)r\Yt)]dt + aldWt 

with given initial perceived emissions E^. Consequently, the aggregate perceived emission E defined 
by 

N 
i=l 

satisfies 

dEt = {bt- f{Yt))dt + atdWt, 

where 

N N N 

h = Y^ bl at = Y, <y\ and fix) = ^[{d)rHx). 

1=1 i=l i=l 

Again, since we are in a Brownian filtration, it follows from the martingale representation theorem 
that there exists a progressively measurable process Z = (^t)tg[o,T] such that 

dYt = ZtdWt and / \Zt\^dt < oo. 

Jo 

Furthermore, in order to entertain a concrete existence and uniqueness result, we assume that W 
is one-dimensional and that there exist deterministic continuous functions [0, T] x M 9 [t, e) ^ 

b(t,e) e R and [0,r] X M 9 t <7(t,e) G M such that bt = b{t,Et) and at = cF{t,Et), for all 
t E [0,r],Q-a.s. 

Consequently, the processes E, Y, and Z satisfy a system of Forward Backward Stochastic Differen- 
tial Equations (FBSDEs for short) under Q, which we restate for the sake of later reference: 

f dEt = {bit, Et) - f{Yt))dt + a{t, Et)dWt, with given EoeR 
dYt = ZtdWt, Yt = X1 [A,+oo) 

The fact that the terminal condition for Yr is given by an indicator function results from the equilib- 
rium analysis of these markets. See |H1 and [H. A is the global emission target set by the regulator 
for the entire economy. It represents the cap part of the cap-and-trade scheme. A is the penalty that 
firms have to pay for each emission unit not covered by the redemption of an allowance. Currently, 
this penalty is 100 euros in the European Union Emission Trading Scheme (EU ETS). Notice that 
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since all the cost functions c* are strictly convex, / is strictly increasing. We shall make the following 
additional assumptions: 

(13) b{t, e) and a{t, e) are Lipschitz in e uniformly in t, 

(14) there exists an open ball U C M^, U 3 (T, A), such thatinf(t e)g[/pi[o,T]xR ^"^{t^ e) > 0, 

(15) / is Lipschitz continuous (and strictly increasing). 

We denote by the collection of all M-valued progressively measurable processes on [0, T] x R, 
and we introduce the subsets: 

Z£lf;E^[ \Zs\'^ds<oo} and := |y G 11°; sup |y,|2] < ooj. 

3.3. Solving the Singular Equilibrium FBSDE. The purpose of this subsection is to prove exis- 
tence and uniqueness of a solution to FBSDE ( 12 1. 



Theorem 1. If assumptions ( [T3| ) to ( [15] ) hold for a given A G M, then, for any A > 0, FBSDE ( [12] ) 
admits a unique solution (E, Y, Z)£S^x§^x M^. Moreover, for any t G [0, T], Et IS non-increasing 
with respect to A and non-decreasing with respect to A. 



Proof. For any function : M M, we write FBSDE(99) for the FBSDE (12i when the function 
g = Al[yv,+oo) appearing in the terminal condition in the backward component of ( 12 1 is replaced by 

(i) We first prove uniqueness. Let [E, Y, Z) and {& , Y', Z') be two solutions of FBSDE ([12 1. 
Clearly it is sufficient to prove that Y = Y'. Let us set: 

6Et := Et - E't, ^Yt ■= Yt - Y^ bZt := Zt - Z^ 

bit,Et)-b{t,E't) ^ a{t,Et)-cT{t,E't) , 

^* SEt MsEtm^ ■- MsEtm- 

Notice that {Pt)o<t<T and {T,t)o<t<T are bounded processes. By direct calculation, we see that 

d{Bt5Et5Yt) = -Bt6Yt{f{Yt) - f{Y^))dt + BtSEtSZtdWt, 

where 



Since SEq = and 5Et5Yt = {Et - E'j){g{ET) - giE^)) > 0, because g is nondecreasing, this 
implies that 

rT 



Bt6Yt{f{Yt)-f{Y:))dt 



< 0. 



Since Bt > and / is (strictly) increasing, this implies that 6Y = dt (>Si dQ— a.e. and therefore 
Y = Y' by continuity. 

(ii) We next prove existence. Let {g"')n>i be an increasing sequence of smooth non-decreasing 
functions with g"' G [0, A] and such that g" — > g^ = Al(yv,oo)- 

(ii-1) We first prove the existence of a solution when the boundary condition is given by g". For 
every n > 1, the FBSDE(5(") satisfies the assumption of Theorems 5.6 and 7.1 in [15] with 63 = 0, 
/i = /2 = /s = 0, CT2 = 0-3 = 0, 62 < (by ([T5])) and h = (since 5" is non-decreasing) so that 
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Condition (5. 1 1) in IBl holds with A = and F{t, 0) = for any e > 0. By Theorem 7.1 in IBl, the 
FBSDE(5") has a unique solution (^", y", Z") eS'^xS'^x if. Moreover, it holds y^" = E]}), 
< t < T, for some deterministic function ti". In contrast with fT5 ], the function is not a random 
field but a deterministic function since the coefficients of the FBSDE are deterministic. We refer to 
ifTTl for the general construction of n" when the coefficients are deterministic. Since the sequence 
{g^)n>i is increasing we deduce from the comparison principle fT5l Theorem 8.6], which applies 
under the same assumption as 1 15, Theorem 7.1], that, for any t G [0, T], the sequence of functions 
(n"(t, •))„>! is non-decreasing. By |15, Theorem 8.6] again, is non-decreasing in A and non- 
increasing in A. Since g"- is [0, A]-valued and e) = K'^[g'^{EJ^)\E^ = e], we deduce that 
is [0,A]-valued as well. Since the sequence of functions (u")n>i is non-decreasing, we may then 
define: 

u{t,e):= lim t^i"(i,e), tG[0,T], eGM. 

n— >-oo 

Clearly, u is [0, A]-valued and u{t, •) is a non-decreasing function for any t G [0, T]. Moreover, u is 
non-decreasing in A and non-increasing in A. 

By Theorem 6.l-{iii) and Theorem 7.1-(i) in |[T5]| . we know that, for every n > 1, the function 
ti" is Lipschitz continuous with respect to e, uniformly in t G [0, T]. Actually, we claim that, for 
any 6 G (0, T), the function u^{t,-) is Lipschitz continuous in e, uniformly in t G [0,T — 5] and in 
n > 1. The proof follows again from Theorem 6.l-{iii) and Theorem in |fT5ll . To be more 

specific, we need to establish a uniform upper bound for the bounded solutions y to the first ODE in 
(T5\ (3.12)] associated with an arbitrary positive terminal condition yx = h > 0. Namely, for given 
bounded (measurable) functions bi : [0,T] B t ^ bi{t) G M+ and 62 : [0,T] 3 t ^ b2{t) G M+, 
with infjgjo b2{t) > 0, we are seeking an upper bound for any bounded {yt)o<t<T satisfying 

yt = yT + J {bi{s)\ys\ -b2{s)y'^s)ds ; yT = h>0. 

Here 61 (t) is understood as an upper bound for the derivative of b with respect to x, and 62 as a lower 
bound for the derivative of / with respect to y. As long as yt doesn't vanish, we deduce from a simple 
computation that 

yt = exp^j^ bi{s)ds^(^^ + &2(s)exp^^ bi{r)dr^ds^ . 

Since the right-hand side above is always (strictly) positive, we conclude that it is indeed a solution 
for any t G [0, T]. Therefore, there exists a constant C, independent of yT, such that yt < C /{T — t) 
for any t G [0, T). By {Hi) in Theorem 6.1 and {i) in Theorem 7.1 in [15], we deduce that, for any 
5 G (0, T], the function •) is Lipschitz continuous with respect to e, uniformly in t G [0, T — 5] 
and n > 1. Letting n tend to +00, we deduce that the same holds for u. 

Notice that the process E"^ solves the (forward) stochastic differential equation 

dE^ = {b{t, E^) - f o u^{t, K))dt + a{t, K)dWt, t G [0, T), 

where here and in the following, we use the notation / o u for the composition of the functions / and 
u. Since / is increasing and the sequence (n")„>i is non-decreasing, it follows from the comparison 
theorem for (forward) stochastic differential equations that the sequence of processes (-E'")„>i is 
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non-increasing. We may then define: 

Et ■■= lim i^t" fort G [0,T]. 



ra— >oo 



(ii-2) To identify the dynamics of the limiting process E, we introduce the process E defined on [0, T) 
as the unique strong solution of the stochastic differential equation 

dEt = {b- fo u){t, Et)dt + cr(t, Et)dWt, t G [0, T) ; Eq = 0. 

The fact that the function u is bounded and Lipschitz-continuous in space (locally in time), together 
with our assumptions on b, f and a guarantee the existence and uniqueness of such a strong solution. 
Since b is at most of linear growth and u is bounded, the solution cannot explode as t tend to T, 
so that the process {Et)o<t<T can be extended by continuity to the closed interval [0, T]. Since u is 
Lipschitz continuous with respect to e, uniformly in t G [0, T — 6] for any 6 G (0, T), we deduce from 
the classical comparison result for stochastic differential equations that E^ > E for any t G [0, T). 
Letting t tend to T, it also holds EJ^ > Et- Since, for any n > 1, (t, e) = [g" {E^)\ = e] , for 
{t, e) G [0, T) X M, and is a non-decreasing function, we deduce that n"(t, .) is a non-decreasing 
function as well. Obviously, the same holds for u{t, •). We then use the fact that > E together 
with the increase of u'^{t, .) to compute, using Ito's formula, that, for any t G [0, T]: 

{K - Et? = 2 - Es) {{b-fo u^)is, E^) -{b-fo u)is, Es))ds 

Jo 

+ I \a{s,E'^)-a{s,Es)\''ds + 2 [ {E^ - Es){a{s, E^) - a{s, Es))dWs 
Jo Jo 

<C [' \E^ -Es\''ds + 2 [\E^-Es){fou-fou^){s,Es)ds 
Jo Jo 

2 [\e^ - Es){a{s,E^) - a{s,Es))dWs 
Jo 

i-t i-t 
<{C + l) \E^-Esfds+ \{f ou- f ou''){s,Es)\^ds 
Jo Jo 



(16) 

+ 



2 [\e^ - Es) {a{s, E^) - a{s, E,))dW,, 
Jo 



by the Lipschitz property of the coefficients b and a. Taking expectation, we deduce 

ft ft 



E'^[{E^ - Etf] < {C + l)E'^ [ \E^-Esfds + E^ [ |(/ o n - / o n")(s, 

Jo Jo 



Then 



E'^[{E'^ - Etf] < + [' E^[{E^' - Esf]ds + e"" 

Jo 

where e" := E'Q [ /(f |(/ o u - / o Es)\ ds] 0, by the dominated convergence theorem. 

Therefore it follows from Gronwall's inequality that supo<^<2- E'^[(£^" — EtY] — )• as n tends to 
+0O. Repeating the argument, but using in addition the Burkholder-Davis-Gundy inequality in ( 16 1, 
we deduce that E"- — > ^ in S^, and as a consequence, E = E. 
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(ii-3) The key point to pass to the hmit in the backward equation is to prove that Q[Et = A] = 0. 
Given a small real 6 > 0,v/e, write 

Q[Et = A] = Q[Et = a, {t,Et)T-s<t<T e U] 
(17) ' ^ -- ^ , 

+ Q[Et = a, 3te[T- 6, T] : {t, Et) U] , 



where U is as in (14i. (Here, the notation {t, Et)T-s<t<T £ U means that {t,Et) G U for any 
t e [T — S,T].) On the event {{t,Et)T-5<t<T G U}, the process {Et)T-5<t<T coincides with 
{Xt)T-5<t<T, solution to 

Xt = ET-5+ f {b{s,Xs)- fou{s,Xs)) + [ a{s,X,)dWs, T-6<t<T, 
JT-5 Jt-5 

where : [0, T] x R R is a given bounded and continuous function which is Lipschitz continuous 
with respect to e, which satisfies inf[o t]xR ct > 0, and which coincides with a on U. 

Since a^^ is bounded and / is bounded on [0, A], we may introduce an equivalent measure Q ~ Q 
under which the process Bt := Wt — d-~^{t, Xt){f o n)(t, Xt), t G [T — 6, T], is a Brownian motion. 
Then X solves the stochastic differential equation 

(18) dXt = b{t, Xt)dt + Xt)dBt, te[T-6,T]; Xt-s = Ets- 

By Theorem 2.3.1 in |[T6l . the conditional law, under Q, of Xt given the initial condition Xts has 
a density with respect to the Lebesgue measure. Consequently, Q[Xt = A] = 0, and the same holds 
true under the equivalent measure Q. Therefore, 

Q[Et = A, {t,Et)T-5<t<T eU]=0. 

By ( fTTj ), we deduce 

?T = A] = Q[Et = A,3tG[T- 6,T] : it,Et) U] 



<q[ sup -(T,^t) I >dist((r, A), ^7' 

T-5<s<T 



As 5 tends to 0, the right-hand side above tends to 0, so that 
(19) Q[Et = A] = 0, 



which implies that we can use g = X1(^a,oo) instead of g = A1[a,oo) in (12 1. Moreover, we also 
have: 

(20) lim Q[^^ > A\Tt] = (Q[Et > MJ't] 

n— >cxD 

for each t < T. The fact that g"" < g implies: 



as n — cxD by (20l. On the other hand, since E'!^ > Et, it follows from the non-decrease of gn, the 
dominated convergence theorem, and (|20j) that 
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Hence Y^" — > Yt := Ef[g{ET)]. Now, let Z € be such that 

Yt = g{ET) - ZsdWs, t E [0, T]. 

Notice that Y takes values in [0, A], and therefore y G S^. Similarly, using the increase and the de- 
crease of the sequences (u")n>i and (i?")„>i respectively, together with the increase of the functions 
u^it, .) and u{t, .) and the continuity of the function u{t, •) for t G [0, T), we see that for t G [0, T): 

u{t,Et) = lim u^{t,Et) < liminf < limsupu"(t, -Ef ) < lim u{t,E^) = u{t,Et). 

n—^oo n—>-co n— >oo n— >oo 

Since Y^^ = u^{t,E^), this shows that Yt = u{t,Et) on [0, T), and the proof of existence of a 
solution is complete. □ 



Impact on the model for emission control. As expected, the previous result implies that the tougher 
the regulation (i.e. the larger A and/or the smaller A), the higher the emission reductions (the lower 
Et)- In particular, in the absence of regulation which corresponds to A = 0, the aggregate level of 
emissions is at its highest. 

We also notice that the assumptions in Theorem[T]can be specified in such a way that the aggregate 
perceived emission process E takes non-negative values, as expected from the rationale of the model. 

Proposition 3. Let the conditions of Theorem^hold true. Assume further that /(O) = and there 
exists r > such that a{t, 0) = 0, b{t, .) > on [0, r], and b{t, .) <0 on [— r, 0]. Then: 

(i) for any Eq > 0, the process E in ( |12[ ) is non-negative; 

(ii) if in addition Eq > 0, then Et > 0/or all t G [0, T). 



Proof. By ( fTS] ), we know that f{y) > for y G [0, A]. Since the process iYt)Q<t<T is [0, A]-valued, 
we deduce from the comparison principle for forward SDEs that the forward process {Et)o<t<T is 
dominated by the solution {Xt)o<t<T to the SDE: 

Xt = Eo+ [ b{s, Xs)ds + [ a{s, Xs)dWs, 0<t<T. 
Jo Jo 

Observe that our conditions on b and a imply that, whenever Eq < 0, we have X < and therefore 
^ < 0. Then Yt = A1[a,+oo)(-^t) = 0, so that n(0, Eq) = E{Yt) = 0. Similarly, n(t, e) = 0, for 
any t G [0, T] and e < 0. ' 

As a consequence, for any initial condition Eq, we can write (/(lt))o<t<T in the forward equation 
in ([T2]l as 

f{Yt) = f{u{t,Et)) = mt,Et)) - f{u{t,0)) = ^^^^''^^^^^"^^^"^''°^^ ^a|^,^o}. 

where the ratio {f{u{t, e)) — f{u{t, 0)))/e, for e 7^ 0, is uniformly bounded in e G M \ {0} and in 
t in compact subsets of [0, T) since u is Lipschitz-continuous in space, uniformly in time in compact 
subsets of [0, T), see Point (ii-1) in the proof of Theorem[T] Similarly, the processes 

_ bit,Et) , , ^ _ a{t,Et) ^ 
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are adapted and bounded, by the Lipschitz property of the coefficients b,a in e uniformly in t, and 
the fact that b{t, 0) = a{t, 0) = 0. We then deduce that {Et)o<t<T may be expressed as 

Et = Eo exp(^ iPs -^s- l^l)ds + ^ ^sdWs^ , 0<t<T, 

with = [fiYt)/Et]l^^^^Qy 0<t<T. □ 

Remark 4. Using for u additional estimates from the theory of partial differential equations, we may 
also prove that ipt appearing in the above proof of Proposition^grows up at most as (T — 
when t T. This implies that ip is integrable on the whole [0, T] and thus, that Ex > Q as well when 
Eq > 0. Since this result is not needed in this paper, we do not provide a detailed argument. 



Remark 5. The non-degeneracy of a in the neighborhood of (T, A), see ( 14 1, is compatible with the 
condition (T(t, 0) = of Proposition^ since A, which is the regulatory emission cap in practice, is 
expected to be (strictly) positive. 



4. Enlightening Example of a Singular FBSDE 



We saw in the previous section that the terminal condition of the backward equation can be a 
discontinuous function of the terminal value of the forward component without threatening existence 
or uniqueness of a solution to the FBSDE when the forward dynamics are non-degenerate in the 
neighborhood of the singularity of the terminal condition. In this section, we show that this is not the 
case when the forward dynamics are degenerate, even if they are hypoelliptic and the solution of the 
forward equation has a density before maturity. We explained in the introduction why this seemingly 
pathological mathematical property should not come as a surprise in the context of equilibrium models 
for cap-and-trade schemes. 



Motivated by the second model given in subsection 2.2 we consider the FBSDE: 



(21) 



with the terminal condition 



dPt = dWt, 

dEt = {Pt-Yt)dt, 

dYt = ZtdWt, < t < T, 



(22) 



Yt 



4A,oo 



for some real number A. Here, {Wt)t^[o^T] is a one-dimensional Wiener process. This unrealistic 
model corresponds to quadratic costs of production, and choosing appropriate units for the penalty 
A and the emission rate e to be 1. (For notational convenience, the martingale measure is denoted 
by P instead of Q as in Section [3| and the associated Brownian motion by {Wt)o<t<T instead of 

iWt)o<t<T)- 

Below, we won't discuss the sign of the emission process E as we did in Proposition [3] above for 
the first model. Our interest in the example (pT])-(22 1 is the outcome of its mathematical analysis, not 
its realism! We prove the following unexpected result. 
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Theorem 2. Given (p, e) G M?, there exists a unique progressively measurable triple {Pt,Et, Yt)Q<,t<T 
satisfying ( |21| ) together with the initial conditions Pq = p and Eq = e, and 

(23) l(A,oo)(£^T)<i^T<l[A,oo)(^T). 

Moreover, the marginal distribution of Et is absolutely continuous with respect to the Lebesgue m.ea- 
surefor any < t < T, but has a Dirac mass at A when t = T. In other words: 

F{Et = A} > 0. 

In particular, (Pt, Et, Yt)o<,t<T may not satisfy the terminal condition F{Yt = 1[a,oo)(^t)} = 1- 
However, the weaker form ( [23] ) of terminal condition is sufficient to guarantee uniqueness. 

Before we engage in the technicalities of the proof we notice that the transformation 

(24) {Pt, Et)o<t<T ^ {Et = Et + {T- t)Pt)o<t<T 
maps the original FBSDE ( [2T] ) into the simpler one 

f dEt = -Ytdt + {T-t)dWt, 
\ dYt = ZtdWt, 

with the same terminal condition Yt = 1[a,oo)(^t)- Moreover, the dynamics of {Et)o<t<T can 



be recovered from those of {Et)o<t<T since {Pt)o<t<T in (21 1 is purely autonomous. In particular, 
except for the proof of the absolute continuity of Et for t < T, we restrict our analysis to the proof 
of Theorem [2] for E solution of (25 1 since E and E have the same terminal values at time T. 



We emphasize that system (25 I is doubly singular at maturity time T: the diffusion coefficient of 
the forward equation vanishes as t tends to T and the boundary condition of the backward equation 
is discontinuous at A. Together, both singularities make the emission process accumulate a non-zero 
mass at A at time T. This phenomenon must be seen as a stochastic residual of the shock wave 
observed in the inviscid Burgers equation 

(26) dtv{t,e) -v{t,e)dev{t,e) = 0, te[0,T),e£R, 



withi;(T, e) = 1[a,+oo)(6) as boundary condition. As explained below, equation ( 26 1 is the first-order 



version of the second-order equation associated with (25 1. 



Indeed, it is well-known that the characteristics of (26 1 may meet at time T and at point A. By 



analogy, the trajectories of the forward process in ( 25 1 may hit A at time T with a non-zero probability. 



then producing a Dirac mass. In other words, the shock phenomenon behaves like a trap into which 
the process {Et)o<t<T (or equivalently the process {Et)o<t<T) may fall with a non-zero probability. 
It is then well-understood that the noise plugged into the forward process {Et)o<t<T may help it to 
escape the trap. For example, we saw in Section[3]that the emission process did not see the trap when 
it was strongly elliptic in the neighborhood of the singularity. In the current framework, the diffusion 
coefficient vanishes in a linear way as time tends to maturity: it decays too fast to prevent almost 
every realization of the process from falling into the trap. 



As before, we prove existence of a solution to ( 25 1 by first smoothing the singularity in the ter 



minal condition, solving the problem for a smooth terminal condition, and obtaining a solution to 
the original problem by a limiting argument. However, in order to prove the existence of a limit, we 
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will use PDE a priori estimates and compactness arguments instead of comparison and monotonicity 
arguments. We call mollified equation the system p5] ) with a terminal condition 

(27) Yt = HEt), 

given by a Lipschitz non-decreasing function (p from M to [0, 1] which we view as an approximation 
of the indicator function appearing in the terminal condition (p2]|. 



4.1. Lipschitz Regularity in Space. 

Proposition 4. Assume that the terminal condition in (|25|) is given by (|27[) w/f/z a Lipschitz non- 



decreasing function (j) with values in [0, 1]. Then, for each (to, e) G [0, T] x M, (25 1 admits a unique 



solution {E^°''^ , yI'^'^ , Zl°''^)tQ<t<T satisfying E^'^'^ = e and Y^'^ = (j)(E!^'^) . Moreover, the map- 



ping 

{t,e)^v{t,e) = Yl'' 

is [0, l]-valued, is of class C^'^ on [0, T) x M and has Holder continuous first-order derivative in time 
and first and second-order derivatives in space. 

Finally, the Holder norms of v, deV, d1 and dtv on a given compact subset of [0, T) x M Jo not 
depend upon the smoothness of (j) provided (p is [0, l\-valued and non-decreasing. Specifically, the 
first-order derivative in space satisfies 

(28) ^<dev{t,e)<^—^, tE[0,r). 

In particular, e ^ v{t, e) is non-decreasing for any t € [0, T). 

Finally, for a given initial condition (to, e), the processes {Yl^^''^)t^-,<t<T ond {Zl'^'^)tQ<t<T> solu- 



tion to the backward equation in ( 25 1 ( with ip as boundary condition ), are given by: 

(29) y/°'' = v{t, ^f''"), to<t<T; = (T - t)dev{t, ^*«'^), to<t<T. 

Proof. The problem is to solve the system 

f dEt = -Ytdt + {T - t)dWt, 
I dYt = ZtdWt, 



(30) 



with ^ = (P{Et) as terminal condition and {to, e) as initial condition. The drift in the first equation, 
i.e. (t, y) e [0, T] X M —y, is decreasing in y, and Lipschitz continuous, uniformly in t. By 
Theorem 2.2 in Peng and Wu |[T8l (with G = 1, /3i = and /32 = 1 therein), we know that equation 



(30 1 admits at most one solution. Unfortunately, Theorem 2.6 in Peng and Wu (see also Remark 2.8 
therein) does not apply to prove existence directly. 

To prove existence, we use a variation of the induction method in Delarue fT|. In the whole 
argument, to stands for the generic initial time at which the process E starts. The proof consists in 
extending the local solvability property of Lipschitz forward-backward SDEs as the distance T — to 
increases, so that the value of to will vary in the proof. Recall indeed from Theorem 1.1 in f7\ 
that existence and uniqueness hold in small time. Specifically, we can find some small positive real 



number 6, possibly depending on the Lipschitz constant of (p, such that (30 1 admits a unique solution 



when to belongs to the interval [T — 6,T]. Remember that the initial condition is Etg = e. As a 
consequence, we can define the value function v : [T — 5,T] x B (to, e) ^ Y-to'^- Corollary 
1.5 in IiTjI, it is known to be Lipschitz in space uniformly in time as long as the initial time parameter 
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to remains in [T — 6,T]. The diffusion coefficient T — tin ( |30l ) being uniformly bounded away from 
on the interval [0, T — 6], by Theorem 2.6 in [7], (30l admits a unique solution on [to, T — 6] when 
to is assumed to be in [0, T — 6). Therefore, we can construct a solution to (30 1 in two steps when 
to < T — 6: we first solve ( 30 1 on [to, T — 5] with Et^ = e as initial condition and v{T — 5, •) 
as giving the terminal condition, the solution being denoted by {Et, Yt, Zt)tQ<t<T-5\ then, we solve 
( 30 1 on [T — 5, T\ with the previous Et-& as initial condition and with (f) as giving the terminal 



condition, the solution being denoted by {E[, y/, Z[)T-5<t<T- We already know that E'j,_^ matches 
Et-5- To patch {Et, Yt, Zt)to<t<T-5 and Y"/, ^t)T-5<t<r into a single solution over the whole 
time interval [to, T], it is sufficient to check the continuity property Yrs 



Y' 



as done in Delarue 



171 . This continuity property is a straightforward consequence of Corollary 1.5 in Q: on \T — 5,T], 



iYl)T-5<t<T has the form Yj. = v{t, E[). In particular, Y,!p_^ - 
Yt-5. This proves the existence of a solution to (30 1 with 



v{T-6,E!j 



T-51 



v{T-6,Et- 



e as initial condition. 

tQ,e -frto,e ryto,e\ 



We conclude that, for any (to, e), (30l admits a unique solution {E°''^, Y^ Z^°'^)tg<t<T satisfy 



y*"'"' (i.e. the 



ing Ef.°' = e and Yj?' = (t){E,^' ). In particular, the value function v : (to,e) ^to 
value at time to of the solution {Yt)tQ<t<T under the initial condition Et^ = e) can be defined on the 
whole [0, r] X M. 

From Corollary 1.5 in Q and the discussion above, we know that the mapping e v{t, e) is 
Lipschitz continuous when T — t is less than 5 and that, for any to € [0,r], y^*0'^ has the form 



y 



io,e 



v{t, E°'^) when T — t is less than 5. In particular, on any [0, T — 6'], 5' being less than 5, 
(30 1 may be seen as a uniformly elliptic FBSDE with a Lipschitz boundary condition. By Theorem 
2.1 in Delarue and Guatteri |9] (together with the discussion in Section 8 therein), we deduce that v 
belongs to C°([0,r] x M) n C^'2([0, T) xM),thatt-^ ||(9ev(t, OIU is bounded on the whole [0,T] 
and that t ^ \\d1^v{t, •)||oo is bounded on every compact subset of [0, T)Q Moreover, (29 1 holds. 



By the martingale property of {Y^ °' )tQ<t<T, it is well-seen that w is [0, l]-valued. To prove that it 



is non-decreasing (with respect to e), we follow the proof of Theorem|lJ We notice that {E^ 
satisfies the SDE: 



to,e\ 



)to<t<T 



dE, 



to,e 



-v{t, El°'^)dt +{T- t)dWt, to<t<T, 



which has a Lipschitz drift with respect to the space variable. In particular, for e < e', e!^'^ < e!^'^ , 



so that v{tQ, e) 



E(t>{E^r°''') < Ecl){E'r° 



v{to,e'). 



We now establish (ESI). For t^ < t < T, the forward equation in (30 1 has the form 



(31) 



El 



v{s,El°''')ds+ / {T-s)dW, 



to 



to 



Since v is in space on [0,T) x M with bounded Lipschitz first-order derivative, we can apply 
standard results on the differentiability of stochastic flows (see for example Kunita's monograph ifTTll ). 
We deduce that, for almost every realization of the randomness and for any t G [to,T), the mapping 



Specifically, Theorem 2.1 in |:9;j says that v belongs to C''([0, T) x E) and that t ^ \\dcv(t, ■)\\oc bounded on every 
compact subset of [0,T). In fact, by Corollary 1.5 in Delarue |7|, we know that v belongs toC°([r - S,T] x R) and that 
t ^ \\dev{t, ■)\\oo is bounded on [T — 5, T] for 5 small enough. 
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eI°'^ is differentiable and 



(32) 

In particular, 
(33) 



deEl"'^ = 1 - f deV{s,El^''')deEl^'^ds. 
Jto 



deEl"'^ = exp 



dev{s,El!''^)ds 



to 



Since v is non-decreasing, we know that deV > on [0, T) x M so that deEl°'^ belongs to [0, 1]. Since 



deV is also bounded on the whole [0, T) x M, we deduce by differentiating the right-hand side in (31 
with t = T that de&r^'^ ^^i^ts as well and that de&T''' = limf^y deEj"''' G [0, 1]. To complete the 
proof of (28 1, we then notice that for any t € [to, T], 

d[{T - t)Yl°'^ - ^*«''] = (T - t)dYl°'^ - (T - t)dWt = {T- t)[Zl'''^ - l]dWt, 

so that taking expectations we get: 

{T - toHto,e) - e = -E[e!^'']. 

Now, differentiating with respect to e, we have: 

(T - to)dev{to, e) = 1 - E[9e^r°'1 ^ 1' 
which concludes the proof of ( [28| ). 

It now remains to investigate the Holder norms (both in time and space) of v, deV, dl^v and dtv. 
We first deal with v itself. For < t < s < T, 

v{s, e) - v{t, e) = v{s, e) - v{s, El'"") + v{s, El''') - v{t, e) 



v{s,e)-v{s,Ein + Y:''-Yi 
v{s, e) - v{s, El'^) + ^ Z^'^dBr. 



From ( [28] ), we deduce 

\v{s, e) — v{t, e)\ < 



< 



< 



1 


T 


— s 




1 


T 


— s 




1 



EE. 



t,e 



e +E 



Zi'^dBr 



T-s 



t + ^ (r - rfdr 
(T - rfdr 



t + 



1/2- 
1/2- 



+ E 



+ {s-t 



iZi'^l'^dr 



,1/2 



1/2 



since Z*'^ = (T — r)dev{r, Ep^) G [0, 1]. So for e > 0, v is 1/2-Holder continuous in time t £ 
[0, T — e], uniformly in space and in the smoothness of <p. 

Now, by Theorem 2. 1 in Delarue and Guatteri ||9l, we know that v satisfies the PDF 



(34) 



dtv{t,e) + 



(T - 



dl^v{t, e) - v{t, e)dev{t, e) = 0, t£ [0, T), e G M, 



with (p as boundary condition. On [0, T — e] x M, e > 0, equation (34i is a non-degenerate second- 
order PDF of dimension 1 with —v as drift, this drift being C^/^'^-continuous independently of the 
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smoothness of (p. By well-known results in PDEs (so called Schauder estimates, see for example 
Theorem 8.11.1 in Krylov HI]), for any small r/ > 0, the C(^-'')/^'^~''-norm of u on [0,T-e] x M is 
independent of the smoothness of (p. □ 



Remark 6. As announced, equation (34i is of Burgers type. In particular, it has the same first-order 
part as equation (|26[). 



4.2. Boundary Behavior. Still in the framework of a terminal condition given by a smooth (i.e. 
non-decreasing Lipschitz) function with values in [0, 1], we investigate the shape of the solution as t 
approaches T. 

Proposition 5. Assume that there exists some real such that (f){e) = 1 on [A^, +oo). Then, there 
exists a universal constant c > such that for any 6 > 

(35) v{t,A+ + T -t + 6) > 1 -exp(-c ^^_^^3 ), < t < T. 

In particular, v{t, e) ^ 1 as t T uniformly in e in compact subsets of{A'^, +00). 

Similarly, assume that there exists an interval (—00, A^] such that 0(e) = on (—00, A~]. Then, 
for any 6 > 0, 

(36) v{t, A^-5)< exp(-c^^--^) . 

In particular, v{t, e) ^ as t T uniformly in e in compact subsets of (—00, A~). 



Proof. We only prove (35 1, the proof of (36 1 being similar. To do so, we fix (to, e) G [0, T) x M and 
consider the following system 

J dE^ = -dt + (T - t)dWt 
\ dYf = Z^rdWt, to<t<T, 

with E^^ = e as initial condition for the forward equation and Yj^ = 4>{E^) as terminal condition 
for the backward part. The solution {eI°'^, Y^"'^, Zl°'^)tf^<t<T given by Proposition |4] with eI°'^ = e 
and y^°'^ = (piE^r^''') satisfies y/"'" € [0, 1] for any t € [to, T] so that E^ < E^^'" almost surely for 
t G [to, T]. Now, since (j) is non-decreasing, (t){E^) < (/)(£'^'^) almost surely, namely Y^ < Y^^"'^ . 
Setting w~(to, e) = Y^^^, recall that Y^^^ is deterministic, we see that: 

(37) v-{tQ,e)<v{tQ,e)<l. 
Now, since 

^;-(to, e) = E</)(^^) = E0^e - (T - to) + - s)dWs^ 
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with (j) > l[A+,+oo)> by choosing e = A+ + {T — to) + 6 as in the statement of Proposition |5] we get: 

E(j){E^) = E(l)(^A+ + 5 + {T-s)dWs^ 

>PA+ + (5+/ (T- s)dWs > A+ 
Jto 

= F I (T - s)dWs >-5 = 1 - P / (T - s)dWs < -6 

and we complete the proof by applying standard estimates for the decay of the cumulative distribution 
function of a Gaussian random variable. Note indeed that var{J^{T — s)dWs) = (T — to)^/3 if we 
use the notation var(^) for the variance of a random variable ^. □ 



The following corollary elucidates the boundary behavior between A and A"*" + (T — t) with A 
and A+ as above. 

Corollary 1. Choose cp as in Proposition^ If there exists an interval [A"*", +oo) on which (j){e) = 1, 
then for a > and e < A+ + {T - t) + {T - we have: 

e - A+ 
v{t,e) > -exp( 



(38) ^..^v 

for the same c as in the statement of Proposition^ 

Similarly, if there exists an interval (— oo, A~] on which (j){e) = 0, then for a > and e > 
A- - (r-t)i+" we have: 



(39) 



e — A , 
v{t,e) < +exp( 



(T-t) 



l-2a. 



+ {T-ty 



Proof. We first prove ( 38 1. Since v{t, ■) is 1/{T — t) Lipschitz continuous, we have: 

A+ - e + (T-t) + (T-f)i+" 



[t, A+ + (T - t) + (T - - v{t, e) < 



A+ 



T-t 



T-t 

+ i + (r-t)". 



Therefore, 



v{t, e) > v{t, A+ + {T-t) + {T- -I- (T-t) 



A+ -e 
T-t 



and applying ( [35| ) 



v{t, e) > - exp(-c(T - tf"-') - (T - ty 



T-t 

For the upper bound, we use the same strategy. We start from 

e-A- 
T-t 



,{t,e)-v{t,A--{T-t)'+") 



< 



+ iT-tr, 
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SO that 

v{t, e) < + exp(-c(r - + (T - t)". 

□ 



4.3. Existence of a Solution. We now establish the existence of a solution to (25 1 with the original 
terminal condition. We use a compactness argument giving the existence of a value function for the 
problem. 



[o,r) X 



[0, 1] satisfying 



Proposition 6. There exists a continuous function v 

(1) V belongs to C^'2([0, T) x M) and solves 

(2) v{t, •) is non-decreasing and 1/(T — t)-Lipschitz continuous for any t G [0, T), 

(3) V satisfies (35 1 and {36 i with = = A, 

(4) V satisfies (38 1 and (39 1 with A" = A+ = A, 

and for any initial condition (to, e) G [0, T) x M, f/je strong solution {Ef^^'^)tQ<t<T of 



(40) 



to 



u(s, Es)ds + {T- s)dWs, to<t<T, 



to 



is such that {v{t, E^°'^))tQ<t<T ci martingale with respect to the filtration generated by W. 

Proof. Choose a sequence of [0, l]-valued smooth non-decreasing functions (i?!>")n>i such that 0"(e) = 
for e < A — 1/n and (/>"(e) = 1 for e > A + 1/n, n > 1, and denote by {v^)n>i the corresponding 
sequence of functions given by Proposition]?] By Proposition [4| we can extract a subsequence, which 
we will still index by n, converging uniformly on compact subsets of [0, T) x M. We denote by v 
such a limit. Clearly, v satisfies (1) in the statement of Proposition 6| Moreover, it also satisfies (2) 
because of Proposition]?] (3) by Proposition [5j and (4) by Corollary 1] Having Lipschitz coefficients, 
the stochastic differential equation (40 1 has a unique strong solution on [to,T) for any initial con- 



dition E, 



to 



e. If we denote the solution by {E'f'^)tQ<t<T, Ito's formula and (34 1, imply that the 



process {v{t, E°' ))to<t<T is a local martingale. Since it is bounded, it is a bona fide martingale. □ 
We finally obtain the desired solution to the FBSDE in the sense of Theorem [2] 
Proposition 7. v and {E'f'^)tQ<t<T being as above and setting 



y/"'" = v{t, El°'^), zl"'" = {T- t)dev{t, ^, 



7^to,e 



to,e 



~?tofi\ 



to<t<T, 



the process {E\^'^)tQ<t<T has an a.s. limit E!^'*^ as t tends to T. Similarly, the process {Y^°''^)t(,<t<T 



?to,<^ 



to,e\ 



has an a.s. limit Yrp°' as t tends to T, and the extended process {Y^ °' )to<t<T <^ martingale with 
respect to the filtration generated by W. Morever, ^-a.s., we have: 



(41) 
and 

(42) 



L(A,oo 



){E 



<y^'"^<i[A,oo)(4°'')- 



y 



to,e 



\rto,'i I 



Z^dW,. 



to 



Notice that Zl°'^ is not defined for t 



T. 
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Proof. The proof is straightforward now that we have collected all the necessary ingredients. We start 
with the extension of {El°'^)tQ<t<T up to time T. The only problem is to extend the drift part in (40 1, 
but since v is non-negative and bounded, it is clear that the process 

/ to<t<T 

is almost-surely increasing in t, so that the limit exists. The extension of (l^/"''^)to<i<r up to time T 
follows from the almost-sure convergence theorem for positive martingales 
To prove ( [4l| ), we apply (3) in the statement of Proposition [5] If e!^ 

't^' ^ ^ - . 



to 



limt-,T Ej"'^ > A, 
I T, 

v{t, El°''') > 1 - exp[-C(5V(r - t)^] for t close to T, i.e. y^°'^ > 1. Since Y^''"" < 1, we deduce 
that 



then we can find some 6 > such that El"'" > A + {T - t) + 6 for t close to T, so that Y^^'" 



E'r^'' > A 



Yr, 



tn,e 



1. 



In the same way, 



E, 



to,e 
T 



< A 



y^°'' = 0. 



This proves (41 1. Finally (42i follows from Ito's formula. Indeed, by Ito's formula and (34l, 



Yto,e ^ yI^'^ + r Zl'^'^dWs, to<t<T. 

Jto 



By definition, Zl°'^ = {T — s)dev{s, Ei^'^), to < s <T.^y part (2) in the statement of Proposition 
[6j it is in [0, 1]. Therefore, the Ito integral 

cT 



to 



makes sense as an element of L {ft, P). This proves (42 1. 



□ 



4.4. Improved Gradient Estimates. Using again standard results on the differentiability of sto- 
chastic flows (see again Kunita's monograph ifTTI ) we see that formulae ( [32] ) and ( [33] ) still hold in 
the present situation of a discontinuous terminal condition. We also prove a representation for the 
gradient of v of Malliavin-Bismut type. 

Proposition 8. For to G [0, T), dev{to, e) admits the representation 



(43) 



WxnviT -5,E^rS^\) 



T 



deEl'^'^'dWt 



to 



In particular, there exists some constant A > such that 
(44) 



sup sup dev{t,e) < +00. 

|e|>ylO<t<T 



Proof. For 6 > 0, Proposition [Vjyields 



E 



{T-6,E'r^fs) / deEl°'^dWi 



to 



zl^'^dWt 



to 
rT-S 



deEl°'^dWt 

I" \T-t)dev{t,El'"')deEl'"'dt 
Jto 
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The bounds we have on deV and {deE^'^' )to<t<T justify the exchange of the expectation and integral 
signs. We obtain: 



E 



v(t-6,e: 



T~S) 



to 



T-S 



to 



{T-t)E[de[v{t,El"^')]]dt. 



Similarly, we can exchange the expectation and the partial derivative so that 



E 



S,E'^'L"s) / deEl''''dWt 

J to 



T-5 



to 



{T-t)de[Ev{t,E;''^^)]dt. 



Since {v{t, E^°''^))tg<t<T-S is a martingale, we deduce 

?to,e 



E 



v(T 



deEl^'^'dWt 



to 



dev{to,e) 
1 



T-S 



(T - t)dt 



to 



{T-5-to){T + 6-to)devito,e) 



Letting 6 tend to zero and applying dominated convergence, we complete the proof of the representa- 
tion formula of the gradient. 

To derive the bound R4h, we emphasize that, for e away from A (say for example e <^ A), the 



probability that {Ef^^'^)tQ<t<T hits A is very small and decays exponentially fast as T — to tends to 



0. On the complement, i.e. for sup^ <j<2^ Ef'^ < A, we know that v{t,E°' ) tends to as t tends 
to T. Specifically, following the proof of Proposition [s] there exists a universal constant c' > such 
that for any e < A - 1 and Iq G [0, T) 

[T - tofdMto, e) < 2{T - to)^/¥i/2 [ sup > A] 

to<t<T 



< 2{T - tof/^F^/^ [A - 1 + sup / (T - s)dWs > A] 

to<t<T Jto 

< 2{T - tof/^F^/^ [ sup / (T - s)dWs > l] 



< 2{T - to)^/'^ ex.p 



the last line following from maximal inequality (IV.37.12) in Rogers and Williams ll20ll . 

The same argument holds for e > A + 2 by noting that (43 1 also holds for t; — 1. □ 

Remark 7. The stochastic integral in the Malliavin-Bismut formula ( |43| is at most of order (T — 
toY^"^. Therefore, the typical resulting bound for dev{t,e) in the neighborhood of (T, A) is (T — 
Obviously, it is less accurate than the bound given by Propositions^and^ This says that the 
Lipschitz smoothing of the singularity of the boundary condition obtained in Propositions ^and^ 
namely dev(t, e) < (T — follows from the first-order Burgers structure of the PDE (34) and that 
the diffusion term plays no role in it. This is a clue to understand why the diffusion process E feels 
the trap made by the boundary condition. On the opposite, the typical bound for dev{t, e) we would 
obtain in the uniformly elliptic case by applying a Malliavin-Bismut formula (see Exercice 2.3.5 in 
Nualart filEl ) is of order (T — t)~^l'^, which is much better than [T — t)~^. 
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Nevertheless, the following proposition shows that the diffusion term permits to improve the bound 
obtained in Propositions^and^ Because of the noise plugged into E, the bound {T — t)^^ cannot 
be achieved. This makes a real difference with the inviscid Burgers equation (|26|) which admits 



(t,e) E [0,r) X 



as solution, with V'(e) = 1 A for e G M. (See for example (10.12') in Lax II13L ) 

We thus prove the following stronger version of Propositions |4] and [6} 
Proposition 9. For any [to, e) G [0, T) x M, it holds (T — to)dev{tQ, e) < 1. 

Proof Given (to,e) G [0,r) x M, we consider (^eI°''' ,¥1°'" , Zl°'%^,<t<T as in the statement of 
Proposition [7] As in the proof of Proposition |4j we start from 

d[{T - t)Yl'''' - El''''] ={T- t)dYl'"' - (T - t)dWt = {T- t)[zl°'' - l]dWt, to<t<T. 
Therefore, for any initial condition (to, e), 

{T-to)v{to,e)-e = -E[E!f>'']. 
Unfortunately, we do not know whether e!^'' is differentiable with respect to e. However, 



{T -to)dev{to,e) = 1 - lim h-'^E[E^ 



to,e+h 



E, 



1 - lim h-^ lim E[£;*°''=+'' - < 1 - lim lim inf E[deE, 



to,e+u~\ 



Using ([33]), the non-negativity of deV and Fatou's lemma. 



(T - to)dev(to, e) < 1 - lim lim inf E 

h^Ot/'T\u\<h 



exp 



h^Ot/'T \u\<h 



to 



< 1 - lim inf E 

/i-s>0 \u\<h 



exp 



to 



< 1 -E 



expf-limsup / 5ev(s,£;*0'"+")d, 

V ^1-^^ \u\<h J to 



Consequently, in order to prove that (T — to)dev{to, e) < 1, it is enough to prove that: 

pto,e+u\ 



(45) 



lim sup / dev{t,Ef''^'')dt 

\u\<hJto 



is finite with non-zero probability. To do so, the Lipschitz bound given by Proposition [4 is not 
sufficient since the integral of the bound is divergent. To overcome this difficulty, we use (|44i: with 



non-zero probability, the values of the process {Et)to<t<T at the neighborhood of T may be made as 
large as desired. Precisely, for A as in Proposition [8] it is sufficient to prove that there exists 6 > 
small enough such that P[inf|/j|<]^ inix-s<t<T 

^to,nh > ^] > 0. For 5 > 0, we deduce from the 

boundedness of the drift in (HOb that 



'[ inf inf E, 

\h\<lT-5<t<T 



to,e+h 



>A\> 



e - 1 - (T - to) + inf 

T-S<t<T 



(T - s)dWs > A 



to 
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By independence of the increments of the Wiener integral, we get 



'hnf inf Ej^'^-^^yA] 

\h\<lT-S<t<T 



> 



e - 1 - (T - to) + 



T-S 



(T - s)dWs > 2 A 



to 



inf 

T-S<t<T 



T-5 



(T - s)dWs > -A 



The first probability in the above right-hand side is clearly positive for T — 5 > to- The second one 
is equal to 

inf / (T - s)dWs > -A =1-F sup / (T - s)dWs > A 

-S<t<T Jrp_g J lT-S<t<T Jt-5 



Using maximal inequality (IV.37. 12) in Rogers and Williams EOl . the above right hand-side is always 

□ 



positive. By ( [44] ), we deduce that, with non-zero probability, the Umsup in ( [45] ) is finite. 
4.5. Distribution of Et for to < t < T. We finally claim: 

Proposition 10. Keep the notation of Propositions^and^and choose some starting point (to, e) G 
[0, T) X M and some p G M. Then, for every t £ [tQ,T), the law of the variable 



obtained by transformation ( 24 1, is absolutely continuous with respect to the Lebesgue measure. At 
time t = T, it has a Dirac mass at A. 

Proof. Obviously, we can assume p = 0, so that Pt = Wt. (For simplicity, we will write E'*"'^ for 



^<o,e,p -J We start with the absolute continuity of El°'^ at time t < T. Since v is smooth away from 



to,e 



T, we can compute the Malliavin derivative of E'*"'^. (See Theorem 2.2. 1 in Nualart lfT6l .) It satisfies 

DsEI°'^ = t-s- I dev{r,El^''^ + {T-r)Wr)DsEl''^^dr- j {T-r)dev{r,E^^'>''' + {T-r)Wr)dr, 

J s J s 

for to < s < t. In particular, 

DsEI"'^ = [1-{T- r)dev[r, + (T - r)Wr)] 



(46) 



u)Wu)du 



dr. 



By Proposition|9| we deduce that DsE^"'^ > for any to < s < t. By Theorem 2.1.3 in Nualart |[T6ll . 



we deduce that the law of i?^*"'^ has a density with respect to the Lebesgue measure. 

To prove the existence of a point mass at time T, it is enough to focus on e!^'"^ since the latter is 
equal to e!^'^. We prove the desired result by comparing the stochastic dynamics of e!^'^ to the time 
evolution of solutions of simpler stochastic differential equations. With the notation used so far, E^^'^ 
is a solution of the stochastic differential equation: 

(47) dEt = -v{t, Et)dt +{T- t)dWt 

so it is natural to compare the solution of this equation to solutions of stochastic differential equations 
with comparable drifts. Following Remark|6j we are going to do so by comparing v with the solution 
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of the inviscid Burgers equation (26 1. To this effect we use once more the function ip defined by 
7/^(6) = 1 A introduced earlier. As said in Remark|7] the function '4>{{e — A)/{T — t)) is a solution 
of the Burgers equation (26 1 which, up to the diffusion term (which decreases to like (T — t)^ 
when t T), is the same as the partial differential equation satisfied by v. Using ( [35] ) and ( [36] ) with 
A" = A+ = A and 5 = {T - t)^/^, we infer that v{t, e) and V(e - A/{T - t)) are exponentially 
close as T - t tends to when e < -(T - tf/^ or e>T-t + {T- tf/^\ using ^ and ^ with 
A~ = A"*" = A and a = 1/4, we conclude that the distance between v{t, e) and ^'((e — A)/(r — t)) is 
at most of order 5/4 with respect to T-t as T-t tends to when -(T-f)^/^ < e < T-t+{T-t)^/^. 
In any case, we have 

A^ 



(48) 



Ve G 



\v{t,e) - 



< C{T - 1 



,1/4 



for some universal constant C. We now compare ( 47 1 with 



(49) 



with X± 



dX~T 



T-t 

e as initial conditions. Clearly, 



)dt ± C(T - tfl^dt + {T- t)dWt, to<t<T, 



(50) 



to<t<T. 



Knowing that ip{x) = x when < 2; < 1, we anticipate that scenarios satisfying < X^ — A < T—t 
can be viewed as solving the stochastic differential equations: 



dZf 



± 



A 



dt ± C{T - tf/^dt + {T- t)dWt, 



with = e as initial conditions. This remark is useful because these equations have explicit 
solutions: 



(51) Z± = A + (T - 1) [Wt - Wt, T 4C7(r - tf'^ ± AC{T - tof'^ + ^ 
We define the event F by: 



A. 



to' 



to<t<T. 



sup \Wt 

to<t<T 



Wto\< 



and we introduce the quantities e{to) and e(to) defined by 



so that 



whenever e(to 



e{tQ 



and 



e(to) 



1 e- 
- < — 
4 - T 



A 3 

— < - 
to " 4 



< e < e(to)- For such a choice of e, since 
Z± - A 



A + -{T-to) 



t 



Wt - Wto ^ AC{T - tf/^ ± AC{T - toY^^ + 



a/4 



e- A 
T-to'' 
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it is easy to see that if we choose to such that T — to is small enough for 32C(T — to)^^^ < 1 to hold, 
then 

on the event F. This implies that {Xf)tf,<t<T and {Zf)tn<t<T coincide on F, and consequently that 
= = A and hence e!^'^ = A on F by (dOI). This completes the proof for these particular 
choices of to and e. In fact, the result holds for any e and any to G [0, T). Indeed, since £^*°'^ has a 
strictly positive density at any time t E (to, T), if we choose ti € (to, so that 32C{T — ti)^/^ < 1, 
then using the Markov property we get 

F{E^r°'^ = A} > / F{E^^'^ = A}P{^*°'' G de'} > 

which completes the proof in the general case. □ 



Remark 8. We emphasize that the expression for DsEj"'^ given in (46) can vanish with a non-zero 
probability when replacing t by T. Indeed, the integral 

cT 



[ dev{u,Ei°'' + {T-u)Wu)du 

J r 



may explode with a non-zero probability since the derivative dev{u,e) is expected to behave like 
(T — u)~^ as u tends to T and e to A. Indeed, v is known to behave like the solution of the Burgers 
equation when close to the boundary, see (48 1. As a consequence, we expect dgV to behave like the 
gradient of the solution of the Burgers equation. The latter is singular in the neighborhood of the final 
discontinuity and explodes like (T — u)~^ in the cone formed by the characteristics of the equation. 

However, in the uniformly elliptic case, the integral above is always bounded since dev{u, •) is at 
most of order (T — u)^^^"^ as explained in RemarkvA 



4.6. Uniqueness. Our proof of uniqueness is based on a couple of comparison lemmas. 

Lemma 1. Let (p be a non-decreasing smooth function with values in [0, 1] greater than 1[a,+oo)' '3«<i 
w be the solution of the PDF ( 34 1 with (p as terminal condition. Then, any solution {E[, Z[)tQ<t<T 
of (25 I starting from E[^ = e and satisfying 1(a,+oo)(^t) — ^t — -'-[A,+oo)(-^t) '^^^^ satisfies 

w{t,E[)>Yl, to<t<T. 
Similarly, if (p is less than 1(a,+oo)> t^^f^ 

w{t,E[)<Yl, to<t<T. 
Proof. Applying Ito's formula to {w{t, E[)tQ<t<T, we obtain 

d[w{t,E[)-Yl] = {w{t,E[)-Yl)dew{t,E[)dt+ [{T - t)dew{t,E'^) - Z[]dWt. 
Therefore, 



[w(t, E'^) - y/] exp(- [ dew{s, E'^)ds 

= exp^- ^ dew{s, E',)ds^ [{T - t)dew{t, E[) - Z[] dWt. 



SINGULAR FBSDES AND EMISSIONS DERIVATIVES 



29 



In particular, 

w{t,E't) -Y^ = expi I dew{s,E'^)ds]E 



to 

which completes the proof. 

The next lemma can be viewed as a form of conservation law. 



□ 



Lemma 2. Let (x")n>i be a non-increasing sequence of non-decreasing smooth functions match- 
ing on some intervals (— oo, A~'")n>i and 1 on some intervals +oo)„>i and converging 
towards 1[a,-|-oo)> ^b,en the associated solutions {w^)n>i, given by Proposition^converge towards v 
constructed in Proposition^ 

The conclusion remains true if{x^)n>i ci non-decreasing sequence converging towards 1(a,+oo)- 



Proof. Each is a solution of the conservative partial differential equation (34i. Considering v"^ as 
in the proof of Proposition |6] we have for any n,m> 1 



{w^-v^){t,e)de= / (x"-,/.-)(e)de, te[0,r). 



Notice that the integrals are well-defined because of Proposition [5] Since (f>"^{e) 
m — )• +00 for e 7^ A, we deduce that 



l[A,+oo)(e) as 



{w'' -v){t,e)de 



4A,+oo) 



;e)]de, tG[0,r). 



Since the right hand side converges towards as n tends to +oo, so does the left hand side, but since 
e) > v{t, e) by Lemma [l] (choosing {E' , Y', Z') = (S*0'^, y*"'^, we must also have: 



lim 

n— >+oo 



|u;"(t,e) - v{t,e)\de = 0. 



Since {w"'{t, •))n>i is equicontinuous (by Proposition |4]l, we conclude that w"'{t, e) — >• v{t, e). The 
proof is similar if l(A,+oo)- C 

To complete the proof of uniqueness, consider a sequence (x")n>i as in the statement of Lemma 
|2] For any solution {E[, Y^, Z[)t^<t<T of (25 1 with E[^ = e. Lemma [l] yields 

w^{t,E't)>Yl, te[to,T). 

Passing to the limit, we conclude that 

v{t,E^)>Y^, te[to,T). 

Choosing a non-decreasing sequence (x")n>i> instead, we obtain the reverse inequality, and hence, 
we conclude that y/ = v{t, E[) for t G [to, T). By uniqueness to (40l, we deduce that E[ = E^"'^, so 



that y/ = y/°' . We easily deduce that Z^ = Zl"'" as well. 

Remark 9. We conjecture that the analysis performed in this section can be extended to more general 
conservation laws than Burgers equation. The Burgers case is the simplest one since the correspond- 
ing forward - backward stochastic differential equation is purely linear 
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5. Option Pricing and Small Abatement Asymptotics 

In this section, we consider the problem of option pricing in the framework of the first equilibrium 
model introduced in this paper. 

5.1. PDE Characterization. Back to the risk neutral dynamics of the (perceived) emissions given 



by ( 12 1, we assume that the emissions of the business as usual scenario are modeled by a geometric 
Brownian motion, so that b{t, e) = be and a{t, e) = ere. As explained in the introduction, this model 
has been used in most of the early reduced form analyses of emissions allowance forward contracts 
and option prices (see [6] and [5] for example). The main thrust of this section is to include the 
impact of the allowance price Y on the dynamics of the cumulative emissions. As we already saw in 
the previous section, this feedback f{Ys) is the source of a nonlinearity in the PDE whose solution 
determines the price of an allowance. Throughout this section, we assume that under the pricing 
measure (martingale spot measure) the cumulative emissions and the price of a forward contract on 
an emission allowance satisfy the forward-backward system: 

Et = Eo+ [\bEs - f{Ys))ds + f aEsdWs 
Jo „ ^0 



(52) 



Yt = \1[a,^){Et)- [ 

K Jt 



T 

ZtdWt, 



with / as in ( [T5] ) with /(O) = and A, A > 0. For notational convenience, the martingale measure is 
denoted by P instead of Q as in Section[3]and the associated Brownian motion by {Wt)o<t<T instead 

of {Wt)o<t<T- 

Theorem [T] directly applies here, so that equation (|52]l is uniquely solvable given the initial con- 



dition Eq. In particular, we know from the proof of Theorem [T] that the solution (lt)o<t<T of the 
backward equation is constructed as a function (Yt = u{t, Et))o<t<T of the solution of the forward 
equation. Moreover, since we are assuming that /(O) = 0, it follows from Proposition [5] that the 
process E takes positive values. 

Referring to lITTl . we notice that the function u is the right candidate for being the viscosity solution 
to the PDE 



(53) 



dtu{t, e) + {be - f{u{t, e)))deu{t, e) + \(T'^e^dl^u{t, e) = 0, (t, e) € [0, T) x M+ 

u(r, .) = Ai[A,oo) 



Having this connection in mind, we consider next the price at time t < r of a European call option 
with maturity r < T and strike K on an allowance forward contract maturing at time T. It is given 
by the expectation 

E{(y,*'^ - K)+} = n{u{r^E'rl - K)^}^ 

which can as before, be written as a function U{t, E^'^) of the current value of the cumulative emis- 
sions, where the notation (t, e) in superscript indicates that Et = e. Once the function u is known 
and/or computed, for exactly the same reasons as above, the function U appears as the viscosity 
solution of the linear partial differential equation: 



(54) 



dtU{t, e) + [be - f{u{t, e)))deU{t, e) + ^a^e^d'^^Uit, e) = 0, {t, e) £ [0, r) x M+ 
UiT,.) = iuir,.)-K)+, 
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which, given the knowledge of u, is a linear partial differential equation. Notice that in the case / = 
of infinite abatement costs, except for the fact that the coefficients of the geometric Brownian motion 
were assumed to be time dependent, the above option price is the same as the one derived in ||51. 



5.2. Small Abatement Asymptotics. Examining the PDEs (53 1 and ( [54| ), we see that there are two 
main differences with the classical Black-Scholes framework. First, the underlying contract price is 



determined by the nonlinear PDE (53 1. Second, the option pricing PDE (54i involves the nonlinear 
term f{u{t, e)), while still being linear in terms of the unknown function U. Because the function u 



is determined by the first PDE (53 1, this nonlinearity is inherent to the model, and one cannot simply 
reduce the PDE to the Black-Scholes equation. 

In order to understand the departure of the option prices from those of the Black-Scholes model, 
we introduce a small parmater e > 0, and take the abatement rate to be of the form / = e/o for 
some fixed non-zero increasing continuous function /q. We denote by n'^ and the corresponding 
prices of the allowance forward contract and the option. Here, what we call Black-Scholes model 
corresponds to the case / = 0. Indeed, in this case, both ([53]) and ([54|) reduce to the linear Black- 



Scholes PDE, differing only through their boundary conditions. This model was one of the models 
used in [5 1 for the purpose of pricing options on emission allowances based on price data exhibiting 
no implied volatility smile. 

For e = 0, the nonlinear feedback given by the abatement rate disappears and we easily compute 
that, for e > 0, 



AE 



l[A,oo)(-£'T 



ln[e exp(6(r - t))/A] a^/T^t 



< t < T, 



(55) n°(t,e) 

(56) [/0(t,e) = E[{u'>{T,E^/n-K) 
where E^'^''^ is the geometric Brownian motion: 

(57) dE'^/'^ = E°/'^ [bds + adWs] , s>t, 

used as a proxy for the cumulative emissions in business as usual, with the initial condition eI'^ = e. 
See for example fSj for details and complements. The main technical result of this section is the 
following first order Taylor expansion of the option price. 



Proposition 11. Let f satisfy ( [T5| ) and {t, e) G [0, r) x (0, +oo). Then, as e ^ 0, we have 

U'{t,e) = U^{t,e) 



eE 



T 



0,t,e 



l[A,oo){u'{T,E^/n) / fo{n^{s,E''/n)deu\sVT,E%'f) 



E' 



0,t,e 



ds 



+ o(e), 



where e ^o(e) — > as e 



0. 



Proof. The proof divided into four parts. 

(i) We first prove that the functions and U^, with u*^ = and C/° = on [0, T] x 
[0, r] X M_ respectively, belongs to C^''^{[0, T) x M) and Ci'2([0, r) x M) respectively. 



and 
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By (|55), we know that n° is C^'^ on [0, T) x M;^. Obviously u° = on the whole [0, T] x {0} 
since A > 0. Using the bound 



exp(-— < W-x exp(-— ), X > 0, 

oo y iTT V " z 



we deduce that 



a(r-t)i/2 |ln(eexp(6(r-t))/A)p, 



" - V^^ln(eexp(6(r-t))/A)|^"P( 2^W^) ^' 

for < e ^ 1, t € [0, T). This shows that u^{t, e) decays towards faster than any polynomial. In 
particular deU^{t,Q) = d1^vP{t,Q) = 0. Differentiating ( [55] ) with respect to e, we conclude by the 
same argument that deU^{t, e) and d1^vP{t, e) decay towards faster than any polynomial, so that the 
first and second-order derivatives in space are continuous on [0, T) x M_(_. Obviously, dtvP{t, 0) = 
for any t G [0, T) and, by differentiating ([55]) with respect to t, we can also prove that dtvP is 
continuous on [0,T) x M+. Since u*^ = on [0,T] x W_, we deduce that v? is of class C^'^ on 
[0, T) X M. 

All in all, the computation of the first-order derivatives yields 

d.uM^e) = |ln(eexp(6(T-t))/A) .^(T-t)/2P 



for e > and t G [0, T). The above right-hand side is less than C /{T — t)^/^ for e away from 0, the 
constant C being independent of t. When e is close to 0, 

\d.u\t.e)\ < ^ (Jln(eexp(6(T-0)/A)P 

so that the bound 

(58) \deU^{t,e)\<C{T -t)-^/"^, tG[0,r), eGM, 

is always true. As a by-product, we deduce that u^{t, E^) < C(T — t)^^/^|£'°|, so that 
{u\t,E^,)-K)+ = (nO(r,i?0) -i^)+l^|^o|>(r_.)i/2K/c}- 

In particular, 

C/°(t,e) = E[{u''{T,E^/n - K)+l^|^o.,e|>(^_^),/,^/^^]. 

By the same argument as the one used for u^, we see that C/*^ and its partial derivatives with respect to 
t and e decay towards as e tends to 0, at a faster rate than any polynomial one. In particular, setting 
[/0(t, e) = for t E [0, r] and e E M, we deduce that C/° belongs to C^'^{[0, r) x M). 

(ii) We use the smoothness of and apply Ito's formula to {u^{s, Es'^'^))t<s<T, where E'^'^'^ 
denotes the forward process in (52]), when / = e/o and under the initial condition E^' = e > 0, 
t E [0, r). Using the fact that u° belongs to C^'^([0, T) x R) together with we deduce that, for 
anyt<S <T, 

rS 
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Clearly, u°{S, e) Al[A,+oo)(e) as 5 / T, for e ^ S. Since F[E'^^'^ = A] = 0, see (ii-3) in the 

proof of Theorem [T] we deduce that (5, ^J^'*'^) — > AE1[a,+oo)(^t*''') = u'{t,e) as S /■ T. 
Therefore, 

(59) u°{t,e)-u'{t,e) = €E fo{u'{s, Ef^^))deU°is, Ef'^)ds, 

the right-hand side above making sense because of the bound ( |58] ). By a similar argument, we get 

(60) U\t,e) = E{u'^{T,Ep^'^) - K)^ + eE [ fo{u'{s, Ef'^))deU^is, El^^''')ds. 



Notice that dell^ is bounded since u^{t, •) is Lipschitz-continuous so that the integral above is well- 
defined. By ( [58] ) and ( |59l ), we know that ||u°(t, •) — u'^{t, ■)\\oo < Ce, for a constant C independent 
of e and r. Therefore, 

E{u\t,E'/'^)-K)^ 

+ E[{u'{r,Ep'n-Ky\l^o^r,EP^.^.Kl<Ce}] 



^{uO{T,E'r-'''')>K+Ce} 



fo{u%s,EI''''))deu\s,Ef^'^)ds 



+ eO{F[\u\T,Ep''^)-K\ < Ce]), 
where O(-) stands for the Landau notation. By ( |60l ), we finally get 



fo{u%s,Efn)deu\s,Efnds 



U^{t,e) = U\t,e) + eE 



fo{u^is,El''n)deu''{s,Ef''^)ds 



(61) 



+ eE y /o(n^(s, EI^'n)deU\s, Ef'')ds 
+ €0{F[\u'^{T,Ep^'^)-K\ < Ce]). 
(iii) We now prove that: 

limP|u°(r, E'/'^) -K\<Ce]=0 
for (t, e) G [0, r) x (0, +oo). For any 6 > 0, 

limsupF[\u^{T,E'/'^) - K\ < Ce] < limsupF[\u°{T, E'/'"") - K\ < C5]. 
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By continuity with respect to parameters of solutions of stochastic differential equations, we see that 
TTTs.i.e — ^ ^o,t,e g _j. Q Therefore, by the porte-manteau theorem, it holds for any 6 > 0, 



that E: 



limsupP[|n°(r,^^'*'^) 



K\ < Ce] < P[|nO(r,£;°'*'^) - K\ < C6]. 



On the interval [K — C6, K + C5], with 5 small enough so that K — C6 > 0, the function u^{t, •) is 
continuously differentiable with a non-zero derivative and thus defines a C^-diffeomorphism. More- 
over, since e > 0, the random variable £'{?'*'^ has a smooth density on the interval [K — C6, K + C5]. 
Therefore, the random variable vP{t^ Er'^'^) has a continuous density on the interval [K — C6, K + 
CS\ . We conclude that 

lim ¥\\u^{t, E^'^n -K\< C5] = 0. 
(iv) We now have all the ingredients needed to complete the proof. From ( [6T] ), we have: 

r '■^ 

C/°(t,e) = U'{t,e) + eE 



fo{u%s,Ef^'))deU^is,El'''')ds 



+ eo(e). 



Since, for any s G [0, T), u'^{s, •) converges towards u^{s, •) uniformly as e tends to 0, and since 

P[nO(r,^°'*'^) 
that 

rT 



K] = 0, we deduce from ([SS]) and from Lebesgue dominated convergence theorem 



U"{t,e) = U'{t,e) + eE 



/ fo{u\s,E''/n)deu\s,E''/nds 



/o {u'{s, E'/n) deU\s, E^/nds 



+ eo(e). 



The final result then follows from the identity: 



deU'^it,e) =E 



E 



}^l^K,o.){u\T,E''/''))deu\T,E'', 



0,t,e\ 



< t < T, e > 0, 



which can be derived by differentiation of ( [56] ) and making use of the equality F[u^{t, Er'^''^) 
K] = 0. 



□ 



5.3. Numerical results. In this final subsection we provide the following numerical evidence of the 
accuracy of the small abatement asymptotic formula derived above: 

(1) We compute numerically n*^ with high accuracy, and we then compute values of If^ using 
the values of u'^ so computed. We used an explicit finite difference monotone scheme (see 
for example lUl for details). The left pane of Figure [T] gives a typical sample of results. For 
the sake of illustration we used the abatement function f{x) = x corresponding to quadratic 
costs of abatement. The penalty, cap, emission volatility and emission rate in BAU were 
chosen as A = 1, A = 1.25, a = 0.3 and b = 2A/T where the length of the regulation 
period was T = 1 year. The prices of the allowances u{t,e) and u'^{t,e) were computed 
on a regular grid in the time x log-emission space. The mesh of the time subdivision was 
At = 1/250. The grid of 1001 log-emission was regular, centered around with mesh Ax 
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Figure 1. European call option prices for e = 0, 0.1, 0.2, • • • , 0.9, 1. 

connected to At by the standard stability condition. We considered an option with maturity 
T = 0.25 and strike K = 0.86. We computed u'^{t,e) and U'^{t,e) over this grid for 11 
values ofe, e = 0,0.1,0.2,-- - ,0.9,1, and we plotted the option prices U'^{t, e) against the 
corresponding allowance prices u'^{t, e). The graphs decrease as e increases from to 1. 

We plotted against n*^ in order to show how the option price depends upon the value of 
the underlying allowance. 
(2) We also computed the expectation appearing as the coefficient of e in the first order expansion 
of Proposition [TT] We used a plain Monte Carlo computation of the expectation with N = 
10, 000 sample paths. The right pane of Figure [T] shows the potential of the approximation 
for e = 0.1. The top plot shows the difference between the exact option value and the linear 
approximation given by setting e = and ignoring the feedback effect. Both option values 
were computed by solving the partial differential equations as explained at the beginning of 
the section. The lower plot shows the first order correction as identified in Proposition 1 1 
showing the potential of the approximation. 

References 

[1] G. Barles and P.E. Souganidis. Convergence of approximation schemes for fully nonlinear second order equations. 

Asymptotic Analysis, 4:271-283, 1991. 
[2] R. Belaouar, A. Fahim, and N. Touzi. Optimal Production Policy under Carbon Emission Market. Preprint. 
[3] R. Carmona, F. Fehr, J. Hinz, and A. Porchet. Market designs for emissions trading schemes. SIAM Review, 52:403- 

452, 2010. 

[4] R. Carmona, M. Fehr, and J. Hinz. Optimal stochastic control and carbon price formation. SIAM Journal on Control 

and Optimization, 48(12):2168-2190, 2009. 
[5] R. Carmona and J. Hinz. Risk neutral modeling of emission allowance prices and option valuation. Technical report, 

Princeton University, 2009. 

[6] M Chesney and L. Taschini. The endogenous price dynamics of the emission allowances: An application to C02 

option pricing. Technical report, 2008. 
[7] F. Delarue. On the existence and uniqueness of solutions to FBSDEs in a non-degenerate case. Stochastic Processes 

and Applications, 99:209-286, 2002. 




36 



RENE CARMONA, FRANCOIS DELARUE, GILLES-EDOUARD ESPINOSA, AND NIZAR TOUZI 



[8] F. Delarue. Estimates of the solutions of a system of quasi-linear PDEs. A probabilistic scheme. Seminaire de Proba- 

bilites, XXXVII, 290-332, 2003. 
[9] F. Delarue and G. Guatteri. Weak Existence and Uniqueness for FBSDEs. Stochastic Processes and Applications, 

116:1712-1742,2006. 

[10] N. Krylov. Lectures on Elliptic and Parabolic Equations in Holder Spaces, volume 12 of Graduate Studies in Mathe- 
matics. American Mathematical Society, 1996. 

[11] H. Kunita. Stochastic Flows and Stochastic Differential Equations. Cambridge Studies in Advanced Mathematics. 
Cambridge University Press, 1990. 

[12] O.A. Ladyzenskaja, V.A. Solonnikov and N. N. Ural'ceva, Linear and quasi-linear equations of parabolic type. Amer- 
ican Mathematical Society, 1968. 

[13] P. Lax Hyperbolic Differential Equations. Courant Lecture Notes. American Mathematical Society, 2006. 

[14] J. Ma and J. Yong. Forward-Backward Stochastic Differential Equations and their Applications, volume 1702 of 
Lecture Notes in Mathematics. Springer Verlag. 

[15] J. Ma, Z. Wu, D. Zhang and J. Zhang. On wellposedness of forward-backward SDEs - a unified approach. Preprint. 

[16] D. Nualart. The Malliavin Calculus and Related Topics. Probability and its Applications. Springer Verlag, 1995. 

[17] E. Pardoux and S. Peng. Backward Stochastic Differential Equations and QuasiUnear Parabolic Partial Differential 
Equations. Lecture Notes in CIS, 176:200-217, 1992. 

[18] S. Peng and Z. Wu. Fully Coupled Forward-Backward Stochastic Differential Equations and Applications to Optimal 
Control. SIAM Journal on Control and Optimization, 37:825-843, 1999. 

[19] H. Pham. Continuous -time Stochastic Control and Optimization with Financial Applications. Stochastic Modelling 
and Applied Probability. Springer Verlag, 2009. 

[20] L.C.G. Rogers and D. Williams Diffusions, Markov Processes and Martingales. Volume 2. 2nd Edition. Cambridge 
University Press, 1994. 

[21] J. Seifert, M. Uhrig-Homburg, and M. Wagner. Dynamic behavior of carbon spot prices. Theory and empirical evi- 
dence. Journal of Environmental Economics and Management, 56: 180-194, 2008. 

[22] Z. Wu and M. Xu. Comparison theorems for forward backward sdes. Statistics and Probability Letters, 79:426-435, 
2009. 

ORFE, Bendheim Center for Finance, Princeton University, Princeton, NJ 08544, USA. 
E-mail address: rcarmona@princeton . edu 

Laboratoire J. -A. Dieudonne, Universite de Nice Sophia-Antipolis, Parc Valrose, 06108 Cedex 02, 
Nice, FRANCE 

E-mail address: f rancois . delarueSunice . f r 

Centre de Mathematiques Appliquees, Ecole Polytechnique, UMR CNRS 7641, 91128 Palaiseau 
Cedex, FRANCE 

E-mail address: gilles-edouard . espinosa@polytechnique . org 

Centre de Mathematiques Appliquees, Ecole Polytechnique, UMR CNRS 7641, 91128 Palaiseau 
Cedex, FRANCE 

E-mail address: nizar . touziSpolytechnique . edu 



